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The quantum phase diagram of the Hubbard chain with correlated hopping is accurately determined 
through jumps in tt in the charge and spin Berry phases. The nature of each thermodynamic phase, 
and the existence of charge and spin gaps, is confirmed by calculating correlation functions and 
other fundamental quantities using numerical methods, and symmetry arguments. Remarkably we 
find striking similarities between the stable phases for moderate on-site Coulomb repulsion: spin 
Peierls, spin-density-wave and triplet superconductor, and those measured in (TMTSF^X. 



The search for electronic mechanisms of superconduc- 
tivity and the study of superconducting and Mott phase 
transitions are among the most interesting subjects of 
the physics of strongly correlated systems. In few cases, 
exact results have helped to elucidate the nature of these 
transitions jj] || . In general one has to rely on numerical 
calculations of finite systems for which quantities like the 
Drude weight D c (which should vanish for an insulator 
in the thermodynamic limit 0), or any other correlation 
function, vary smoothly at the transition. Consequently, 
for instance, the boundaries between a charge-density- 
wave (CDW) or spin-density-wave (SDW) insulators and 
metallic phases in half-filled generalized Hubbard models 
were difficult to establish |)]||. 

The Berry phase is a general geometrical concept which 
finds realizations in various physical problems 0. It is 
the anholonomy associated to the parallel transport of a 
vector state in a certain parameter space. In condensed 
matter, the charge Berry phase j c is a measure of the 
macroscopic electric polarization in band or Mott insula- 
tors H while the spin Berry phase 7 S represents its spin 
polarization pjicj|. In systems with inversion symmetry 
7 C and j s can attain only two values: or tt (mod(27r)). 
Thus, if two thermodynamic phases differ in the topolog- 
ical vector 7 = (7c, 7 S ) this sharp difference allows us to 
unambiguously identify the transition point even in finite 
systems. This "order parameter" was recently used to de- 
tect metallic, insulator and metal-insulator transitions in 
one-dimensional lattice fermion models pjlt| . 

In this Letter we determine the quantum phase dia- 
gram of the Hubbard chain with correlated hopping at 
half-filling using topological transitions. The phase di- 
agram is very rich showing two metallic and two insu- 
lating thermodynamic phases each characterized by one 
of the four possible values of the topological vector 7. 
One of the metallic phases corresponds to a Tomonaga- 
Luttingcr liquid with dominant triplet superconducting 
correlations at large distances (TS). This is interesting 
since there is experimental evidence indicating that the 
Bechgaard salts (TMTSF%C1C>4 and (TMTSF)PF 6 un- 
der pressure are TS O-nJ. Furthermore, the insulat- 



ing SDW and spin gapped spin-Peierls phase observed 
in (TMTSF)PF 6 as the pressure is lowered M are also 
present in the model phase diagram. 
The effective model Hamiltonian is: 0] 

H= ^ (4* c j* + h - c -) {pAA (1 - n is )(l - n js ) 

+ tsB n is n ]S + t AB [n is (l - n ]S ) + n js (l - n lS )]} 
+ uY,( n *-\)(nn-\). (1) 

i 

H contains the most general form of hopping term de- 
scribing the low energy physics of a broad class of sys- 
tem Hamiltonians in which four states per effective site 
are retained. In particular, the Hamiltonian H in Eq. 
(1) has been derived and studied for transition met- 
als, organic molecules and compounds 0, intermedi- 
ate valence systems, cuprates and other superconduc- 
tors p5| . In the continuum limit, the only relevant in- 
teractions at half-filling are U and Iaa + ^bb — Hab 
Ifrjf . Therefore, we restrict the present study to the 
electron- hole symmetric case (t^A = Ibb = 1) which 
has spin and pscudospin SU(2) symmetries, the latter 
with generators ij + = Yli(~ l) 8 c]^.c|^, T]~ = (?/ + ) t , and 
rf = ^X^(Ecr n io ~ !)■ The canonical transformation 
(CT) c,-f = Cif, Cii — {—l) l c\^ changes the sign of U in 
H , and interchanges the total spin and pseudospin oper- 
ators {rj a < — ► S a ). These symmetry properties become 
crucial in this work. For tAB =0, the model has been 
solved exactly || with the result that the ground state 
(GS) is highly degenerate. For tAB ^ the physics of the 
model is still unclear and constitutes our main concern. 
For the present case 7 ClS are defined as || 

7c, s = i I # (9k(4>, ±<t>)\d^g K {<t>, ±0)} , (2) 
Jo 

where \gK{<Ph 4>l)) is the GS in the subspace with total 
wave vector K and other quantum numbers kept fixed, 
with fluxes <\> a for spin a. Changes in macroscopic polar- 
ization with spin cr, P a , are related to the corresponding 
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changes in the Berry phase by: APf ± AP± = eA r y CtS /2w 
(mod(e)) Thus, a phase transition will be de- 

tected by a jump in tt of j c ("f s ) if and only if both ther- 
modynamic phases differ in P| + Pj (Pf — P±) by e/2 
(mod(e)). For example, if one of the phases is a CDW 
with maximum order parameter (CDWM) and the other 
a Neel state (N) , one is transformed into the other trans- 
porting half of the charges (those with a given spin) one 
lattice parameter. In addition, as explained below, in 
the present model topological transitions in j c and 7 S 
indicate the opening of the charge and spin gap A c , A s . 

We find that the minimum of the GS energy as a func- 
tion of fluxes, E g (<fii , <fii), corresponds to the so-called 
closed shell conditions (CSC): if the number of sites (as- 
sumed even) is L = 2 (mod(4)), then K = § a = 0, while 
for L multiple of four K = <p a = n (which is equivalent 
to taking antiperiodic boundary conditions and K = in 
a system without fluxes). 
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FIG. 1. Quantum phase diagram of the correlated hopping 
Hubbard chain. The vector Berry phase y = (7 C ,7 S ) and the 
nature of each stable phase is indicated: Luther-Emery liquid 
with equally decaying singlet superconductor and CDW cor- 
relations (SS-CDW), Tomonaga-Luttinger liquid with triplet 
superconducting and bond SDW correlations dominating at 
large distances (TS-BSDW), SDW insulator (SDW), and spin 
and pseudospin dimerization (S,PS). 

For any finite system and fixed tAB, varying U, two 
topological transitions occur in the model, correspond- 
ing to a jump in either 7 C or 7 S . We have determined 
those transitions in rings of length L = 6, 8, 10, 12 using 
the Lanczos method. The results extrapolated with a cu- 
bic polynomial in 1/L are represented in Fig. [jj In con- 
trast to other physical quantities which show large finite- 
size effects, particularly near tAB = |(|, the topolog- 
ical transitions converge rapidly to the thermodynamic 
limit (for example, for tAB = 0.05, "f c jumps at [7=3.451, 
3.681, 3.788, and 3.846 for L = 6 - 12, and the ex- 
trapolated value is [7=3.932). The numerical conver- 
gence becomes problematic for smaller tAB values. At 
tAB = the transition points are determined from the 
exact solution as those values of U where A c and A s 



open. Those critical values are U c>8 — ±4 and match 
smoothly with the rest of the curves in Fig. [l[ It is easy 
to see that under CT the geometrical phases transform 
as 7 C < — ► 7 S + 7r 0. Thus, as seen in Fig. |l|, a jump in 
7 C at U c (full line) implies a jump in j s at —U c (dashed 
line), and vice versa. 

In the case where all particles are localized one can 
easily determine the value of 7 as 7 C , S = Imlnz^ 8 , where 

z l = \9\e ^ J Iff) , (3) 

was recently used to study quantum localization [|l7],|l8| . 
In the thermodynamic limit z c L vanishes for a conductor 
while |z£| — > 1 for an insulator. Clearly, 7(CDWM) = 
(0, 0), while 7(N) = (tt, tt), and by continuity -f(CDW) = 
(0,0), 7 (SDW) = (-7T, 7r). On the other hand, it is not easy 
to predict the values of 7 in conducting phases. However, 
for [7 = 0, the model is invariant under CT. Therefore 
7 = (tt, 0), or 7 = (0, 7r), indicating a topological differ- 
ence with the above mentioned CDW and SDW states. 
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FIG. 2. Spin gap A s as a function of U for tAB = 0.6, 
obtained with DMRG. 

A c (A s ) vanishes at the left (right) of the full (dashed) 
line in Fig. ^, and is different from zero at the right 
(left). It has been shown, that in spin SU(2) symmet- 
ric systems the opening of A s can be detected as a level 
crossing of triplet and singlet states for boundary condi- 
tions opposite to the CSC ones (periodic if the number 
of particles with a given spin N a is even, and antiperi- 
odic if N„ is odd) Il9] . It is precisely this crossing which 
causes the jump in 7 S po[ . A direct evaluation of A s has 
large finite-size effects, and an accurate calculation re- 
quires use of the Density Matrix Renormalization Group 
(DMRG) method Q. Fig. | displays A s for t AB = 0.6 
extrapolated from calculations of open chains with length 
L < 40. For — 2 < U < — 1, A s vanishes within numer- 
ical accuracy (~ 0.01) while for smaller values of [7 it 
increases rapidly. These results are consistent with A s 
opening near U s « —2.1 with a singular growth (expo- 
nentially small as in the t — J model Jl9]|). The result 
derived from the topological transition is U s = —2.051. 



2 




Fig. 3 Aligia 



FIG. 3. Drude weight D c and z c L for tAB = 1.4 as a function 
of U for L — 12 (open circles) and L = 24 (solid triangles). 
The solid lines are polynomial extrapolations in 1/L. 

The symmetry transformation CT implies that if A s 
opens at U s , a pseudospin gap A,, opens at —U s . Pseu- 
dospin excitations form a subset of all charge excitations. 
In our case the charge velocity v c and A c (computed as 
usual in finite systems L < 12 [^9|) coincide with their 
pseudospin counterparts v n and (computed from v s 
and A s for opposite U) [Ell. This is consistent with the 
exact solution for tAB — where the charge excitations 
of lowest energy are pseudospin ones ||. In addition, 
the charge-charge and spin-spin correlation functions are 
interchanged as U changes sign (see below). The open- 
ing of A c where 7 C jumps from to 7r is also consistent 
with calculations of z£, D Cl superconducting correlation 
functions, K p and central charge c. For tAB > 1, A c 
opens more slowly and the detection of the transition 
becomes more difficult. However, as shown in Fig. [|, D c 
and z c L display a similar behavior near the jump in 7 C 
(£/ c =-1.702) as in the quarter-filled infinite U extended 
Hubbard model as a function of the nearest-neighbor re- 
pulsion V p8[ , where a metal- insulator transition takes 
place at V — 2t. At large \U\ the only relevant energy 
scale is 4t\ B /\U\ and therefore, D c increases with U for 
large negative U. Near U c there is a drastic change of 
behavior of D c and zl vs U, and for U > U c , the extrap- 
olated values suggest a tendency to reach the insulating 
values (D c =0, z£=l), in the thermodynamic limit. z£ 
always decreases (increases) with L at the left (right) of 
U c . This is the first accurate (DMRG) calculation of z£ 
in a system of more than 16 sites. 



In order to further characterize each thermodynamic 
phase we use symmetry arguments and the following cor- 
relation functions (CF) (see Fig. 4): 

Xt(s){d) = 2~(( c ot c H c oj. c iT)( c ci+U c <2T ^~ c d+iT c dl)) 
Xos{d) = (c^cjjc^c^) 

Xc(d) = -^(("OT + "-ox - l)(n dT + n dl - 1)) 

Xbsdw{d) = 2(( c ot c U + c iT c oi)( c d+U c dT + c di c d+i-\)) 
Xsz (d) = -2(S* S* d ) , (4) 

computed with DMRG and CSC |l). The CF are defined 
such that, in the non-interacting case, lun^^^ x{d, L) = 
l/(7T<i) 2 , with d an odd number. Xbsdw{d) corresponds 
to correlations between spins located in bonds [fl6|| . 
The GS for CSC is always a spin and pseudospin sin- 
glet. Invariance under pseudospin rotations implies that 
Xos(d) = (— l) d+1 Xc(d), i.e. on-site singlet superconduct- 
ing and charge CF have the same distance dependence. 
This, in turn, implies K p — 1 in the conducting phases 
p2|| , meaning that the large distance behavior is dom- 
inated by logarithmic corrections. Also, using CT and 
spin rotations one obtains Xc(d,U) = Xsz(d, —U) and 
X t(d, U) = X t(d, -U) = (-l) d x bsdw (d, U). 

We restrict our study to odd values of d for which the 
oscillatory factors of the CF in the continuum limit are 
maximum p2[ . In addition we took L — 2d, motivated 
by results on the Heisenberg model p2j ], showing that 
lmif^oo Xsz(d, L) = Cxsz(d,2d), where C is a constant 
independent of d. This also true, if d is odd, for all CF of 
our model in the non-interacting case, with C — (ir/2) 2 . 
The above mentioned symmetry relations and analytical 
results in the non-interacting limit allowed us to check 
the accuracy of the DMRG results. A summary for each 
non-equivalent topological region is described below: 

a) 7 = (tt,0), A e ^ 0, A s ^ 0. All CF decay exponen- 
tially. To understand the nature of this insulating GS let 
us analyze the limit tAB —> oo, where effectively there is 
a sequence of spin and pseudospin dimers. Each dimer is 
the ground state of the model for two sites and two par- 
ticles. Including the dimer-dimer interaction in second- 
order perturbation theory the resulting GS energy per 
site is e w -1. 0625^5 (cdmrg = -1.0808^ B ). The 
energy difference between the two lowest singlet states 
(with K = and K = ir) decays nearly exponentially 
with L, indicating a breaking of the translational sym- 
metry in the thermodynamic limit. 

b) 7 = (ir, 7r), A c 7^ 0, A s = 0. As in the ordinary pos- 
itive U Hubbard model, the GS is an insulating SDW. 

c) 7 = (0,0), A c = 0, A s ^ 0. The system forms a 
Luther-Emery liquid with Xs and Xc CF decaying as 1 jd 
(neglecting logarithmic corrections.) 

d) 7 = (0,7r), A c = 0, A a = 0. Xos and Xs display the 
same long distance behavior (Fig. El). Because of sym- 
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metry arguments, for [7 = 0, all CF except Xt = Xbsdw 
show the same 1 /d 2 decay (apparently without logarith- 
mic corrections). Renormalization group arguments p^ ] 
show that when A s = 0, Xt should decay more slowly 
than Xs, and therefore dominate for small \U\. As U in- 
creases, Xt decays more rapidly, while the opposite hap- 
pens with Xsz- Near the opening of A c (at U c w 2.05 for 
tAB = 0.6), both CF seem to decay in a similar fashion. 
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FIG. 4. Correlation functions (CF) times the square of the 
distance d as a function of d for tAB = 0.6: nearest-neighbor 
triplet \xt(d)\ = \Xbsdw(d)\ (squares) and singlet Xs(d) (cir- 
cles) pair CF, on-site pair and charge CF |x s(d)| = |Xc(°OI 
(downward triangles), and spin CF Xsz(d) (upward triangles). 
See Eq. (§). 

In conclusion, we have constructed the quantum phase 
diagram of Eq. (Q) from topological considerations. Each 
thermodynamic phase is associated to a topological vec- 
tor, and changes in that quantity signal the transition 
point. A key finding is the identification of a triplet 
superconducting phase degenerate with a bond-located 
SDW for tAB < 1 and small \U\. Taking into account 
the slight dimerization in (TMTSF)aX compounds, an 
effective Hamiltonian similar to Eq. (jl]) can be real- 
ized, where only four low-energy states per unit cell are 
kept. Assuming that additional interactions stabilize the 
TS with respect to the BSDW it is remarkable that, for 



U > 0, as tAB is decreased a similar sequence of quantum 
phase transitions takes place as pressure is applied [ jl4| . 
This general approach can be extended to spatial dimen- 
sions higher than one j^], and applied to more general 
models which do not necessarily have SU(2) pseudospin 
symmetry 0. 
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